Spin-Peierls-like transition in AFe 2 As 2 (A=Ba, Sr) 
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From first-principles density functional theory calculations combined with varying temperature 
Raman experiments, we show that AFe2As2 (A=Ba, Sr), the parent compound of the FeAs based 
superconductors of the new structural family, undergoes a spin-Peierls-like phase transition at low 
temperature. The coupling between the phonons and frustrated spins is proved to be the main cause 
of the structural transition from the tetragonal to orthorhombic phase. These results well explain 
the magnetic and structural phase transitions in AFe2As2(A=Ba, Sr) recently observed by neutron 
scattering. 

PACS numbers: 74.25.Jb, 71.18.+y, 74.70.-b, 74.25.Ha 



Introduction The recent discovery of iron-based high- 
transition temperature (high-T c ) superconductors [l| has 
invoked great research interests in similar materials with 
Fe-As layers. Currently, the focuses are mainly on two 
kinds of structures: the first is iron arsenide-oxides with 
P4/nmms space group 0,0,0], represented by the par- 
ent compound LaFeAsO; and the second is ternary iron- 
arsenide compound with body-centered I4/mmm space 
group [H, 0,0] , represented by the parent compound 
BaFe2As2[5(. Those compounds show similar structural 
and magnetic properties. At room temperature, they all 
lay in non-magnetic, high-symmetric state, with Fe-As 
layers separated by La-0 layers or Ba layers respectively. 
With the decrease of temperature, LaFeAsO undergoes 
a slight structural transition from the tetragonal to or- 
thorhombic phase at T a = 150K, followed by the ap- 
pearance of a magnetic SDW(spin-density wave) state 
at 134K0, while BaFe2As2 experiences a similar struc- 
tural and magnetic transition simultaneously at about 
T a = 140K @. At the magnetic SDW state, they all 
form a collinear stripe-ordering magnetic ground-state, 
with the nearest Fe atoms aligning anti-ferromagnetically 
along one crystal axis in Fe-As plane, while ferromagnet- 
ically parallel to the other axis SHE EH- By doping 
with electron or hole carriers, the structure transition and 
the SDW state are both suppressed, and the supercon- 
ductivity emerges at 52K[i|] and 38K[l| respectively. 
Thus there exists a competition between magnetism and 
superconductivity in these compounds, but the subtle de- 
tails such as the pairing mechanism are far from clear. 

As that in high-T c cuprates, revealing the mecha- 
nism of superconductivity of iron-based superconductors 
highly requires understanding the electronic, structural 
and magnetic properties of parent compounds first. The 
structural and magnetic phase transitions seem to be 
common features in these Fe-As based superconductors, 
thus revealing the sources of these transitions and the 
possible connections between them may lead to better 
understandings of the experimental observations. Some 
first-principles density functional theory (DFT) calcula- 



tions have suggested that the frustrated superexchangc 
interactions between Fe ions induces the collinear stripe 
antiferromagnetic ground state [l^, 14 1, on the other 



hand the main cause of the structural phase transition, 
as well as its correlation to the magnetic one, is not well 
understood so far. 

Let us start with the frustrated J\ — J2 Heisenberg 
model to describe the nearest neighbor and next-nearest 
neighbor superexchange interactions among the Fe atoms 
bridged by As atoms, which can be described as [n| 



H 



j\ Si ■ Sj + J2 2_, &i ■ Sj 



(i) 



<%3> 



where < ij > and <C ij S> denote the summation over the 
nearest and the next-nearest neighbors respectively. The 
ground state of the frustrated J\ — J2 spin-half model on 
a square lattice at zero temperature has been studied by 
several groups in the literature and the main results are 
summarized as follows 15j, [l6[: (1) In the absent of frus- 



tration (J2 = 0), its ground state has long-range Neel or- 
der; (2) With the increase of frustration (J2/J1), a phase 
transition from Neel order to a spin-liquid phase occurs; 
(3) If further increasing the frustration, a collinear order 
emerges at J%/ J\ > 0.55, with the nearest spins aligning 
anti-ferromagnetically along one axis while ferromagnet- 
ically parallel to the other. 

According to the DFT calculation in Ref. llfj, J\ almost 
equals to J2 for Fe-As superconductors. In this sense, 
the DFT and model calculations consistently explain the 
collinear magnetic order of the parent compounds. How- 
ever, the ground state given by J\ — J2 model is twofold 
degenerate with the 7r/2 rotational symmetry, which is 
not in agreement with the orthorhombic structure ob- 
served in experiments. It indicates that the pure spin 
model is not sufficient to account for the structural phase 
transition. We hereby suggest to extend J\ — J2 model by 
evolving the spin-phonon coupling. For examples, when 
adiabatic phonons are considered, the model should be 
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modified as, 



H = £ [Ji(l - aiyij)Si ■ Sj 



K 



1 .2 ■ 



£ [Ml-a^S.-Sj + ^y^l (2) 



where 



Uj| with Uj denoting the in- plane dis- 



placement of atom i, a is the spin-phonon coupling con- 
stant. The ground state (at T = 0) of above model was 
calculated with a spin-wave approximation in recent lit- 
erature and a Peierls-like transition from a tetragonal 
to an orthorhombic phase was found at large frustration 
(J 2 /Ji> 0.5)0. 

In this letter, motivated by the model analysis, we use 
the first-principles density functional calculations com- 
bined with varying temperature Raman experiments to 
study the ground state of AFe 2 As 2 (A=Ba, Sr). We 
prove that the spin-Peierls-like phase transition is the 
very mechanism of the structural transition at T a in the 
parent compounds of iron-based superconductors. 

Method The calculation was done using a plane- 
wave based method(l8| with local spin density approx- 
imation (LSDA) and generalized gradient approxima- 
tion (GGA) of Perdew-Burke-Ernzerh (PBE)p]| for the 
exchange-correlation potentials. The density-functional 
perturbation theory (DFPT) was used to calculate 
the T-point phonons. Firstly, the electronic proper- 
ties of AFe2As2(A=Ba,Sr) were calculated to deter- 
mine the electronic ground state, using experimental 
cell parameters [H and energy- minimized internal atomic 
positions. Different magnetic configurations, namely 
the nonmagnetic, square anti-ferromagnetic and collinear 
anti- ferromagnetic were considered |10|. The ground 
states of both BaFe2As2 and SrFe2As2 are found to be 
collinear anti-ferromagnetic with a slight structure tran- 
sition, which is in agreement with previously reported 
results. For the phonon calculations, we focus on the non- 
magnetic high-temperature state and the collinear anti- 
ferromagnetic ground state at low temperature. Here we 
adopt the triclinic primitive unit cells as shown in Fig. [IJ 
in which the structure transition alters the angle between 
x and y crystal axes from 90° to 89.6°. 

The Raman measurements were performed with a 
triple-grating monochromator (Jobin Yvon T64000), 
which works with a microscopic Raman configuration. A 
50 x objective microscopic lens with a working distance 
of 10.6 mm, is used to focus the incident light on sample 
and collect the scattered light from sample. The detector 
is a back-illuminated CCD cooled by liquid nitrogen. An 
solid-state laser (Laser Quantum Torus 532) with high- 
stability and very narrow width of laser line, is used with 
an excitation wavelength of 532 nm. The laser beam of 
3 mW was focused into a spot of less than 10 microns in 
diameter on sample surface. 




FIG. 1: (color online) Structure of BaFe2As2, from left to 
right: conventional cell in I4/mmm space group, primitive 
non-magnetic cell and primitive stripe-ordering cell 



Result From symmetry analysis, the Raman phonon of 
nonmagnetic AFe2As2 with I4/mmm space group con- 
sists of four modes: the Ai s and E s for As, the Bi g and 
Eg for Fe. The vibrating directions in A± g and Bi ff modes 
are perpendicular to the Fe-As layer, while that of E g par- 
allel to the Fe-As layer. After cooling down to the tran- 
sition temperature T a , the crystal structure changes into 
orthorhombic F/mmm space group, and the magnetism 
appears almost at the same temperature. In this space 
group, the E s modes split into B2 ff and B3 9 modes, and 
Ai s mode for As changes into B s mode. The calculated 
Raman modes, as listed in Table 1, are in good agree- 
ment with symmetry analysis. All of the modes have 
been found and assigned from atomic displacements. 

The Raman modes of SrFe2As2 have been measured 
in Reffiol but only for nonmagnetic state reported. By 
comparing our calculated phonon frequencies with the 
experimental values, we find a systematic frequency shift 
(~ 20cm" 1 ) of the calculated values towards higher 
value region (see Table 1) due to the temperature ef- 
fect. That is to say, the DFPT theory deals with the 
ground-state problems, which in this case corresponds to 
a non-magnetic state at T = OK, while the Raman spec- 
tra was measured at room temperature, as a consequence, 
the corresponding experimental values are lower. 

Since the strong spin-phonon interaction is the basis of 
the suggested mechanism of spin-Peierls-like phase tran- 
sition, we must verify this point first, by calculating the 
change of the Raman modes in nonmagnetic state with- 
out the magnetic transition but with the structural tran- 
sition. In that case, the calculated Raman modes and 
frequencies almost remain unchanged for both BaFe2As2 
and SrFe2As2, with frequency shifts (or splits) less than 
3 cm -1 . This suggests that the slight structure transi- 
tion has limited effects on the electronic properties and 
Raman phonons. Our tentative calculations also gives us 
another valuable information: the E g mode of Fe atom 
experiences a frequency splitting as a result of the struc- 
tural transition, which consists with the fact that the 
structure distortion mainly appears in Fe-As layer and 
changes the nearest neighbor Fe-Fe distances slightly, so 
that the nearest neighbor Fe-Fe distances are no longer 
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equal in the orthorhombic crystal axes directions. 



Table 1: Raman phonon frequencies of AFe2As2(A=Ba,Sr) in the non-magnetic state (marked as N) and in the 
collinear anti-ferromagnetic state (marked as C). The Raman modes are in I4/mmm space group corresponding to 
the non- magnetic state, while those in the brackets are in P/mmm space group for the collinear anti- ferromagnetic 

state. The measurement temperatures are 260K and 87K for SrFe2As2, and 290K and 100K for BaFe2As2- The 
atomic displacements arc for the collinear anti-ferromagnetic state. 



Atom 


Raman Mode 


SrFe 2 As 2 -N 


SrFe 2 As 2 -C 


BaFe 2 As 2 -N 


BaFe 2 As 2 -C 


Displacement 




(cm -1 ) 


Cal. (Exp.) 


Cal. (Exp.) 


Cal. (Exp.) 


Cal. (Exp.) 


of Atoms, C 


As 


E s( B 2 S ) 


138.9 


96.1 


140.2 


86.8 


As(x), Fe(y) 


As 


E 9 (B 3g ) 


138.9 


126.3 


140.2 


125.6 


As(y), Fe(x) 


As 


Ai 9 (B 9 ) 


207.6 


183.4 


205.7 


177.8 


As(z) 


Fe 


Blg(Blg) 


219.5(206) 


219.3(213) 


224.0(206) 


214.2(212) 


Fe(z) 


Fe 


E ff (B 2 g) 


301.2 


262.4 


293.6 


252.0 


Fe(y), As(x) 


Fe 


E s( B 3g) 


301.2 


289.0 


293.6 


281.1 


Fe(x), As(y) 



We are now on the position to elucidate the changes 
of the Raman modes due to spin-phonon interactions 
at the magnetic ground state. Our calculations on 
the electronic structure verify that the magnetic ground 
state of AFe2As2(A=Ba,Sr) is a collinear stripe-ordering 
anti-ferromagnetic one. The Fe spins align parallel 
along the shorter axis in Fe-As layer and anti-parallel 
along the longer one, as observed by neutron diffraction 
measurement [Sj], which is similar to that of LaFeAsOQ. 
To correctly illustrate this magnetic configuration, the 
nonmagnetic primitive unit cell should be doubled in a- 
b plane to include two formula cells. Here we adopt a 
-\/2x-\/2xl cell, with x(y) axis rotates 45° and points to 
the nearest Fe-Fe directions. By doubling the crystal lat- 
tice, the reciprocal lattice rotates and shrinks to a half, 
and the original reciprocal points fold down to the new 
points. The folding is illustrated in Fig. [2] (the structure 
distortion from tetragonal to orthorhombic not shown in 
the figure). It is clear that original M(tt, it, 0) points 
fold down to the new F(0, 0, 0) point. Thus the obtained 
r point phonon modes in SDW state doubly contain the 
information from both original T and M points. Only the 
modes corresponding to T point are picked up according 
to the atomic displacements, which are shown in Table 1 
(marked as C). For the purpose of comparison, we also 
show the measured Raman phonons at different temper- 
atures before and after the magnetic transitions in the 
table. The calculated phonon frequencies are found very 
close to the experimental values, which further proves the 
reliance of our phonon calculations. 

As shown in Table 1, accompanied with the onset of 
magnetism, almost all of the Raman phonon frequencies 
softened to lower values, specifically some of them (As- 
B2 S and Fe-B2 ff ) show a giant phonon softening (GPS). 
The GPS of Raman phonon modes is the main result of 
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FIG. 2: Reciprocal space of BaFe2As2 in non-magnetic(solid 
lines) and stripe-ordering(dashed lines) state. 

our DFT calculations, which directly support the mech- 
anism of the spin-Peierls-like phase transition at T a in 
AFe2As2(A=Ba, Sr). Let us elaborate this point in more 
details as follows: 

Firstly, by comparing the changes of B 3g and B 2g fre- 
quencies of Fe atoms (seen in Table 1), one can see that 
the former only reduces 12cm _1 for both BaFe2As2 and 
SrFe2As2 while the latter reduces 41cm _1 for BaFe2As2 
and 39cm~ 1 for SrFe2As2 (GPS). The difference between 
these two modes is the involved Fe atoms vibrating in 
different direction - along x and y crystal axis respec- 
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tively. In SDW state, the spins on Fe atoms align an- 
tiparallel along x (longer) direction, while parallel along 
y (shorter) direction. In B 2 (, mode, the vibration direc- 
tion of Fe atoms accords with the parallel spin direction, 
suggesting an instability of parallel alignment of Fe spins. 

Secondly, if assuming y<y> = 5 x ' y , = S x + 

8 V , ot\ — ai — a and K\ = Ki = K in the extend 
Ji — J2 model [see Eq. (j^J], one can obtain the collincar 
state (ground state) with S x — a(J 2 — Ji)/4K and S y — 
a(J 2 + Ji)/AK at sufficient large frustration Jij J\ > 0.5 
17| . In consideration of J2 ~ J\ in AFe2As2, one can 
reach S v 3> 5 X , which consistently explains why the GPS 
of Raman modes of Fe atoms mainly happens along y 
crystal axis. 

Thirdly, the spin-Peierls-like phase transition proved 
by the DFT and model analysis here closely relates 
to the "exchange striction" effect that also predicts 
symmetry breaking distortions from the magnetoelastic 
interaction 2l| . That effect is believed to play an impor- 
tant role for the origin of the multiferroics, and further- 
more, we prove here it may be also related to the high-T c 
superconductivity. 

Summary In summary, the consistence of the first- 
principles density functional calculations and model anal- 
ysis strongly prove that the spin-Peierls-like phase tran- 
sition is the very mechanism of the structural transition 
at T a in parent compounds of iron-based superconduc- 
tors. We thus conclude that the dominate interactions 
in those parent compounds are frustrated spin-spin in- 
teraction and spin-phonon interaction. Electron or hole 
doping suppresses the ground state of parent compounds 
(orthorhombic structural and collinear magnetic phase) 
and induces the superconductor state. In order to under- 
stand the pairing mechanism for ion-based new supercon- 
ductors, one should treat the spin-spin and spin-phonon 
interactions on an equal footing. 

Supports from the National Natural Science Foun- 
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National Basic Research Program of China (Grants 
No. 2007CB925001) are gratefully acknowledged. 
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